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ABSTRACT 

Analytical formulae for triple differential distributions d 3 cr/ (dcos#dsids2) in the 
neutral current process e + e~ — > /1/1/2/2 are given. They allow to obtain angu- 
lar distributions, rapidity distributions and transversal momentum distributions 
of one fermion pair by only two numerical integrations. Cuts can be applied to 
the integration variables. 



1. Introduction 



Many events of four fermion final states are expected to be observed at LEPII 
and NLC (next linear collider). The accuracy of the theoretical description of such 
processes should match the future experimental precision. As is known from LEPI, 
both, semianalytical programs and Monte Carlo generators are needed to understand 
the physics of experimental data. 

Semianalytical calculations lead to fast FORTRAN codes with high numerical pre- 
cision. Although a semianalytical calculation will never be able to simulate a real 
detector, it should be as differential as reasonable. This allows to predict observables 
closer to the experiment and to make more meaningful comparisons with other groups 

The semianalytic approach has been applied earlier to the calculation of four 
fermion final states of neutral current process i e + e~ — > /1/1/2/2 and Higgs produc- 
tion in the same reaction @. There, analytical formulae for the double differential cross 
section were given. The total cross section was obtained by two numerical integra- 
tions. Differential distributions over the invariant energy of one fermion pair fi) 
were obtained by one integration. 

In this contribution, we present analytical formulae for triple differential cross 
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sections. In addition to the invariant energies of the two fermion pairs, we keep the 
cosine of the angle between the 3-vector of the momentum of one fermion pair and 
the beam axis as an additional parameter. The resulting analytical formulae turn 
out to be not much longer than the double differential ones. However, the formulae 
presented here may be used to get new distributions by two numerical integrations: 
angular distributions , rapidity distributions and transversal momentum distributions 
of one fermion pair. These distributions can be obtained including additional cuts on 
the remaining integration variables. 

We start with the process-independent description of the phase space in section 2 
and give formulae for the cross section and distributions in section 3. Section 4 
contains a brief summary. Detailed formulae for the kinematical functions of the 
triple differential distributions are given in the appendix. Numerical applications will 
be presented elsewhere. 

2. Phase space 

We are interested in the cross section and in distributions of the reaction 

e + (h)e-(k 2 ) -> /i(pi)/i(p 2 )/ 2 (P3)/2(p4) (1) 

withp 2 = m 2 , i = l,...,4. 

The eightdimensional phase space of the four particle final state is parametrized 

as 

i=l A Pi z=l 

yJX(s,s u s 2 ) j\{si,m\,ml) y/\(s 2 ,ml,ml) 

= 2n- dsids 2 di«ii2idiZ2- (2) 

8s 8si 8s 2 

A is the kinematical function, 

A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2ac - 26c, A = A(s, s±, s 2 ), (3) 

and the invariants s, si and s 2 are 

s = (h + k 2 ) 2 , si = (pi + p 2 ) 2 , s 2 = (p 3 + Pif- (4) 

We integrate over every spherical angle fl (fli, fl 2 ) of the e + e~ pair (j\f\, f 2 f 2 
pair) in its rest frame. The spherical angles are decomposed in two integrations, 
dfli = dcos9id(fii with the kinematical ranges 

-l<cos0i<l, < <pi < 2tt. (5) 
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In particular, we have df2 = dcos^d^, where cos 6 = c is the cosine of the angle 

— # 

between the 3- vectors (pi + p 2 ) and k\. In the case without transversal beam po- 
larization, the integration over the rotation angle around the beam axis is trivial, 
giving 2n. 

In the following, we will integrate the squared matrix element over the four re- 
maining angular variables cos 9i, 4>i, cos 9 2 and 2 . This results to analytic functions 
depending on the remaining variables c, s\ and s 2 . They enable us to calculate not 
only total cross sections but also angular distributions da/dc, transversal momentum 
distributions da/dpix or rapidity distributions dcr/dy, of the "compound particle" 
(fi, fi, i = 1,2). These new distributions are useful for comparisons with Monte 
Carlo programs and for more reliable experimental predictions. 

The transformation of the phase space volume to the new integration variables is 

rvtrtv 1 ) rs%(piT,yi) 



f 1 f s t f s t^) _ fVi fPiriV 1 ) fsJ(piT,yi) 

/ dc / _ dsi /_ ds 2 = 2a/s / dy dpi T ds 2 

J-l Js 1 Js 2 (si) Jy 1 Jp 1T (vi) Js 2 (piT,yi) 

fPtr fyfiPir) rs%(piT,Vi) 

= 2Vs / dp 1T dyi / _ ds 2 (6) 

J p 1T J y 1 (piT) Js 2 (piT,yi) 

with similar formulae for p 2 r and y 2 . The old variables c and Si have to be expressed 
through the new variables pyr and y: 



c = 



1 +pj T /( m i + m 2) 2 
\ 1 + coth 2 (y) p\ T l{m-y + m 2 ) 2 ' 



Si = 2 v / sy(m 1 + m 2 ) 2 +p? T coshy + s 2 - s (7) 
The kinematical ranges of the integration variables s\ and s 2 {s\) are 



s 1 = (m 1 + m 2 ) 2 , sf = (y/s - m 3 - m 4 ) 2 , 
sj(si) = (m 3 + m 4 ) 2 , s^(si) = - v^i) 2 - 



(8) 



For simplicity, we give the other kinematical borders only for the case of massless 
final fermions: 

y - = -oo, yt = oo, 

PiAvi) = °> Pit(Vi) = V^/coshy, 

Pit = 0, pt T = y/s, ^ ' 

yf(Pir) = -arcosh^, yf(pi T ) = arcosh^| 

The borders of s 2 (yi,pir) are 
«2 (VuPit) = maxj 0, s - 2 v / ipirC0shy| , s 2 (yi, Pit) = (y/s - Pit cosh y) 2 . (10) 
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Fig. 1: The crab diagrams. 



3. Cross sections and distributions 



We present the formulae for the triple differential distribution of the reaction (|TJ), 
where only neutral gauge bosons are exchanged. The initial electron-positron pair 
and the final fermions f\ and f 2 are all different and none of them belongs to one 
electroweak multiplet. The reaction (|l|) is then described by 24 Feynman diagrams 
which exchange photons and Z-bosons. In the case of four quarks in the final state, 
there are 8 additional diagrams with gluon exchange. One can divide these diagrams 
in three different types. The first type which we call crab (these diagrams are also 
called conversion diagrams) is shown in fig. 1. These diagrams can be considered as 
e + e~ — > fifi with initial state radiation of a /2/2-pair (or e + e _ — > f 2 f 2 with initial 
state radiation of a /1/1-pair). The second type which we call deers (they are also 
called annihilation diagrams) can be considered as e + e~ — > f\fi with radiation of 
a /2/2-pair from the final state (/i-deer), see fig. 2. A different set of diagrams is 
obtained for e + e~ — > f 2 f 2 with radiation of a /1/1-pair from the final state (/ 2 ~deer). 
The process is symmetric in f\ and /a- 

The squared matrix element of the considered reaction is calculated using the 
symbolic manipulation program FORM i. As a result, analytical formulae for the triple 
differential cross section are obtained, 

d 3 a(c, s, s 1 , s 2 ) _ y, d 3 a k (c, s, si, s 2 ) ^ 
dcdsids2 dcdsids2 

The summation runs over the squares and interferences between the three different 
types of diagrams: k — 1 corresponds to the square of crabs, and k = 2, 3 to that 
of /i-deers and / 2 ~deers. Further, k = 4, 5 and 6 correspond to the interferences 
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Fig. 2: The /i-deer diagrams. The / 2 -deers may be obtained by interchanging f\ 
and f 2 . 



between /i-deers and / 2 -deers, between crabs and / 2 -deers, and between crabs 
and /i-deers, respectively. 

For the squares of diagrams of one group (k = 1, 2, 3), we obtain 

d 3 o-i(c, s, si,s 2 ) v 7 ^ , , , , n - rc ( v 

t~] — j — — n^422ie, s; ji, Si; j 2 , s 2 iiy 422 {c; s; s 1 , s 2 , 

dcdsids 2 7rs 

dV 2 (c,s,si,s 2 ) (f , s . psVrC c / x 

"T~ i j — K ( ~'422{Jl,Si,J 2 ,S 2 ,e,Sl 2 y 422 {C,Si,S 2 ,S , 

dcdsids 2 ns 2 



d 3 a 3 (c,s,s 1 ,s 2 ) 
dcds!ds 2 



^C 422 (/ 2 , s 2 ; e, s; /i, si)T 3 ^4 22 (c; s 2 ; s, si) . (12) 



The color factor is T k = N c (fi)N c (f 2 ) for all k = 1, . . . , 6, if photons or Z-bosons are 
exchanged. N c (f) is equal to unity for leptons and three for quarks. In the case of 
two different quarks in the final state, also gluon exchange is possible. The interfer- 
ence between diagrams with one gluon and those without a gluon gives %. — 0, the 
interference between two diagrams with a gluon gives T k = 2/9, k = 2,3, 4. The cor- 
responding contributions have to be added to a 2 , a 3 and <r 4 . The kinematical function 
Q% 22 {c\ s; si, s 2 ) is given in the Appendix. The function C 422 contains couplings and 
gauge boson propagators, 



C422(e,s;/i,si;/ 2 ,s 2 ) 



2s 1 s 2 



3?e 



E 



( 67r2 ) 2 vwfcy^* D vM D vM D vM) D vM*) 
x [L{e,V^L{e,V k )L{e,Vj)L{e,Vi) 

+R(e,V i )R(e,V k )R(e,V j )R(e,V l )} 
x [L(/i, V^)L(/i, V fc ) + V fc )] 
x [L(/ 2 , V-)L(/ 2 , V,) + i2(/ 2 , VS-)i2(/ 2 , VI)] . (13) 



We use the following conventions for the left- and right-handed couplings between 
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vector bosons and a fermion /: 

L(f,l) = R(f,j) = e -^, 

L(q,g) = R{q,g) = -y/4^r B . (14) 
The propagators are 

D v (s) = s-M*+i^r v (s), (15) 

where My and IV(s) = y^Fy/My are the mass and width of the exchanged gauge 
boson (M 7 = M g = T 7 (s) = T g (s) = 0). Note that ai(c, s, s\, s 2 ) is of special interest 
because it describes the production and decay of two off-shell Z-bosons, see fig. 1. 
This is the numerically largest contribution, if the final fermions are not both quark 
pairs. 

For the interferences between different groups of diagrams (k = 4, 5, 6), we obtain 



d 3 a A (c, s, si, s 2 ) 




dcdsids 2 


TCS 2 


d 3 a 5 (c, s, s 1 ,s 2 ) 


= Vx 


dcdsids2 


ITS 2 


rfV 6 (c, s, s 1 ,s 2 ) 


= ^ 


dcdsids 2 


ITS 2 



C 233 (f2, s 2 ; e, s; / x , Si)7^ 33 (c; ^2; s; si) . (16) 

Again, the kinematical function {? 233 (c; s 2 ; s, si) is given in the Appendix. For C 2 33 
we get 

r> ( f t \ 2sSlS2 so V- 1111 

( 67r ) v^ky^z D Vi(s) D Vj (s 2 ) D* Vk (s) D Vi ( Sl ) 

x [L(e,Vi)L(e,V k ) + R(e,Vi)R(e,V k )\ 

x [L(fi, Vi)L{h, V 3 )L(f u Vi) - R(fx, V^Rtfi, ^)i?(/i, V t )] 

x [L(f 2 , V 3 )L(f 2 , V k )L(f 2 , V t ) - R(f 2 , V 3 )R(f 2 , V k )R(f 2 , V t )} . 

(17) 

The triple differential distributions are described by only one combination of cou- 
plings for every interference, exactly as the double differential case i. However, the 
kinematical ^-functions are different. They depend on one more variable c which 
singles out the electron-positron pair and makes the ^-functions less symmetric. 

We are now prepared to write the formulae for various distributions, 

da r s l . /-4( s i) , d 3 a(c, s, Si, s 2 ) 

dsi / ds 2 - 



dc Js-( Sl ) dcdsids 2 



da fPrrM . r4(PiT,yi) d 3 a(c, s, s 1 , s 2 ) 
— = / dp 1T / ds 2 — , 

da rvtiviT) r4(piT,yi) d 3 a(c, s, si, s 2 ) 

1 = / d ^ / dS2 A A A • 

dpiT Jviipvr) J 82 (pit, vi) dcds 1 ds 2 
The kinematical borders of integration and the transformations from c and s% to y\ 
and pit are given in (0) - ([10]) . 

The distributions d 3 a(c, s, si, s 2 )/(dcdsids 2 ) are derived for massless fermions. 
Hence, the formulae are not applicable in the phase space regions si ~ (mi + 
m 2 ) 2 , s 2 ~ (rn.3 + m.4) 2 . However, the mass effects can be completely removed by a 
moderate cut on si and s 2 . Such a cut should be applied to quark pairs in any case 
to remove non-perturbative QCD effects. To leading order, mass effects can be taken 
into account as in ref. @. 

The numerically most important radiative corrections to the considered process 
are due to initial state radiation of photons. They may be taken into account with a 
convolution formula which is well known from the Z line shape !'!. 



4. Summary 

In this paper, compact analytical formulae for triple differential cross sections 
d 3 a(c, s, Si, s 2 )/(dcdsids 2 ) of the process QTJ) are presented. They make the semiana- 
lytic approach more flexible allowing the calculation of many distributions by at most 
two numerical integrations. Further, it allows the inclusion of cuts on the integration 
variables as summarized in the table: 



observable 


cuts possible on 


No. of num. integrations 


a 


Sl, s 2 


2 


da/dsi 


s 2 


1 


da/dc 


Sl, s 2 


2 


da/dpiT 


yi,s 2 


2 


da/ dyi 


Pit, s 2 


2 



Table 1: Examples of cross sections and distributions calculable by the semianalytic 
approach. 

In addition to the observables indicated in the table, double differential distri- 
butions can also be calculated by less numerical integrations. The triple differential 
distribution d 3 a/(dcdsids 2 ) and the double differential distribution d 2 a/(dsids 2 ) can 
be used as input for further investigations. 

The formulae presented here, can easily be generalized to the exchange of extra 
heavy neutral gauge bosons by an extension of the summation in C4 22 and C 2 ^. The 
generalization to longitudinal polarized beams is straightforward. 
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Appendix A Formulae for the kinematical functions 

We start with the kinematical function for crab diagrams squared, see fig. 1: 

„CC, s -1/1 1 \ S 2 + ( Sl + S 2 ) 2 If I 1\ 1 /A „ N 

gg^c;s; Sl ,s 2 ) = — [- + -) _ V y - + 8l8a - (A.l) 



with 



Ti = -\{^ ~ s i - s 2 + cy/X^j , T 2 = -- (s - si - s 2 - cVX V (A. 2) 



The fi~deer diagrams squared lead to the function 

3 

(c;si;s 2 , s) = -(1 + c 2 )Q 42 2{su sz, s) 



+ ^—^-s 1 (s + s 2 )^- (l - 2C( Sl ; s 2 , s) ^ ) . (A.3) 

A 4 V si — s 2 — sj 

The logarithm £(s; si, s 2 ) is defined as 

£(s; si, s 2 ) = —?= In -=. (A.4) 

V A s — Si — s 2 — v A 

^422 (si; S2, s) is the corresponding kinematical function for the double differential 
distribution d ^ 2 ^ which we give here for sake of completeness i: 



g 422 (s; Si, s 2 ) = £(s; Si, s 2 ) - 2 (A. 5) 

S - Si - s 2 
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The function of / 2 -deers squared is given by G^{c\ s 2 ; S\, s) which expresses the 
symmetry of the problem according to the two final fermion pairs. 

We can unify the two functions flA.lj ) and ( |A.3| ) in one function (/£ 22 (c; Si; s 2 , s) 
to be close to the old notation used for double differential distributions: 



£422(0; s;si,s 2 ) 



£422 ( c ! s 'i s ii s 2) if s > si and s > s 2 
(c; s; si, s 2 ) else. 



(A.6) 



As in the double differential case, G% 22 {c\ s i s i> s 2) is symmetric in the last two argu- 
ments. 

The kinematical function of the interference between j\-deers and f 2 -deers is 



£233 (c;s;si,s 2 ) 



;(1 + C 2 )£?233(>; Si, S 2 ) 

3 3 



A 2 8 



(1 -3c 2 )s 



£(si; s 2 , s)2s 2 (si - s 2 ) + (s - s 1 - 3s 2 ) 



x 



£(s 2 ; s, si)2si(s 2 - si) + (s - s 2 - 3si) 



• (A-7) 



£233(5; si, s 2 ) is the kinematical function for the double differential distribution H: 



£233(5; si, s 2 ) 



— j£(s 2 ;s,si)£(si;s 2 ,s) 



4s 



ssi(s - s x ) 2 + ss 2 (s - s 2 ) 2 + SiS 2 (s a - s 2 )' 



+ (s + Si + s 2 ) 



£(s 2 ; s, si)2s (s - s 2 ) 2 + si(s + s 2 - 2s x ) 



+ £(si; s 2 , s)2s (s - si) 2 + s 2 (s + Si — 2s 2 



+ 5s 2 - 4s(si + s 2 ) - (si - s 2 ) ; 



(A.8) 



The interference between crabs and /i-deers is given by: 



1/1 



£233 (c; S2; s; si) = ^S — + )£(si;s 2 ,s 



1 



x4s 2 
+ (s + Si + s 2 ) 



ssi(s - si) 2 + ss 2 (s - s 2 ) 2 + sis 2 (si - s 2 y 



[s - s 2 ) 2 + si(s + s 2 - 2si 



-1/1 1\ 

T k + tJ ^ 

+ £(s x ; s 2 , s)s 2 (si - s 2 ) 2 + s(si + s 2 - 2s) 
+ i[5s 2 -4s 2 (s + si) - (s-si) 2 



(A.9) 



9 



The interference between crabs and f 2 -deers is expressed by ^23? ( c > s i5 s > s z)- ^ n 
contrast to the double differential case, Q 23 ^(c; s ^ s > s i) * s n °t symmetric under the 
exchange of the last two arguments. Again we write the functions ( |A.7| ) and ( |A.9| ) 
formally as one: 

nc (r . „. q . „ ^ _ / G 23 f(c;s;s u s 2 ) if s > si and s > s 2 , . 

C, 233 (c, s, s h s 2 ) - I ^£» (c . g . Si . S2) elge _ (A.lUj 

All kinematical functions are finite and independent of c for A — ► (unfortunately, 
in ref. i is a typing error in the formula for the limit of ^233(5; si; S2)): 

rc , \ s{s 1 + s 2 ) 



/ \ si + s 2 -5s . 

^ 33(C;S;S1;S2) ~> {S1 - S2 - S){s2 - S - Sl) ( A ' U ) 

Finally, we remark that the above functions give those of the double differential 
distributions, if integrated over c: 

J i dc ^4 22 (c; s; s x , s 2 ) = £422(5; s x , s 2 ), 

J i dc G 233 (c; s i *i; s 2) = fe3(s; si, s 2 ). (A.12) 
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